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' Abstract 

We extend our earlier investigation of the finite temperature deconfinement transition in 
SU(iV) gauge theories, with the emphasis on what happens as iV — > oo. We calculate the 
latent heat, L^, in the continuum limit, and find the expected behaviour, Lh oc iV 2 , at large 
N. We confirm that the phase transition, which is second order for SU(2) and weakly first 
^ ■ order for SU(3), becomes robustly first order for N > 4 and strengthens as N increases. As an 
aside, we explain why the SU(2) specific heat shows no sign of any peak as T is varied across 
what is supposedly a second order phase transition. We calculate the effective string tension 
and electric gluon masses at T ~ T c confirming the discontinuous nature of the transition 
for N > 3. We explicitly show that the large- N 'spatial' string tension does not vary with 
T for T < T c and that it is discontinuous at T = T c . For T > T c it increases oc T 2 to a 
good approximation, and the fc-string tension ratios closely satisfy Casimir Scaling. Within 
very small errors, we find a single T c at which all the fc-strings deconfine, i.e. a step-by-step 
breaking of the relevant centre symmetry does not occur. We calculate the interface tension 
but are unable to distinguish between the oc iV or oc iV 2 variations, each of which can lead to a 
striking but different N = oo deconfinement scenario. We remark on the location of the bulk 
phase transition, which bounds the range of our large- iV calculations on the strong coupling 
side, and within whose hysteresis some of our larger- N calculations are performed. 



1 Introduction 



In previous papers, [HE], we have presented results on some basic properties of the SU(iV) 
deconfining phase transition in 3 + 1 dimensions, with the emphasis on the order of the 
transition and the value of the deconfining temperature T c . In the present paper we shall 
address a number of more detailed issues such as the physical latent heat, the mass gap and 
interface tensions. We will also take the opportunity to use some new calculations to improve 
upon our earlier calculations of T c /y/a where a is the confining string tension. In a companion 
paper |3 3 we have provided a detailed study of the properties of topological fluctuations across 
the phase transition. 

In the next Section we provide a background discussion to motivate the calculations that 
appear in subsequent sections. We then briefly summarise our (entirely conventional) lattice 
setup, go on to our improved analysis of T c /y/a and then proceed to our calculation of the 
physical latent heat and its iV-dependence. In the following section we present various mass 
calculations, both at T ~ T c and at higher T. We then address the issue of whether there 
are multiple deconfining transitions at larger N. We follow this with an estimate of the N- 
dependence of the confined-deconfined interface tension. We then return to SU(2) to examine 
more carefully why the usual thermodynamic correlators show no sign of a phase transition 
in that case. Finally we summarise some properties of the bulk phase transition and how it 
affects our calculations. In the concluding section we summarise what we have learned in this 
study, with some comments on large- N Master Fields and Eguchi-Kawai space-time reduction. 

2 Deconfinement 

In this Section we shall briefly summarise some theoretical expectations for the quantities that 
we shall calculate later on in this paper. 

We begin with a very basic discussion of deconfinement and of the dynamical significance 
of the order of the transition. We then discuss a new phenomenon that may arise for N > 3: 
deconfinement may occur through several phase transitions rather than just a single one. 
Although this seems possible on general grounds we shall argue that for dynamical reasons it 
is very unlikely to occur in the case of SU(iV) gauge theories. We then discuss the volume 
and N dependence of the deconfining transition and how this depends on the iV-dependence 
of the latent heat and of the confining/deconfining interface tension. We describe some of 
the strikingly simple scenarios that may result at N = oo. We then discuss our expectations 
for the masses which we calculate below, above and at the phase transition. These include 
effective string tensions, spatial string tensions and the electric Debye mass. Finally we remind 
the reader of the fact that not every aspect of the phase transition on our Euclidean space-time 
volumes is directly reflected in the physics of the hot gauge theory. This will motivate our later 
investigation of the SU(2) phase transition where, as is well known (at least to practitioners), 
there is no peak visible in the specific heat despite the apparent second order nature of the 
transition. 
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2.1 Introduction 



We begin by recalling the following simple argument for deconfmement jlj. Consider two 
static sources in the fundamental representation that are a distance r apart. Suppose that 
the theory is linearly confining, so that the flux is restricted to a flux tube of finite width 
that ends at the two sources. Clearly such a flux tube will be of length I > r. Consider the 
case where / ^> r. The thermodynamic weight of such a configuration depends on the string 
energy E(l) = al, where a is the string tension, as well as on the density of states, i.e. on the 
number, n(l), of flux tubes of length I with fixed end points: 

probability oc n(l)e t 1 oc e d e~^ . (1) 

Here we use the fact that for I ^> r (and ignoring various inessential complications) n(l) ~ 
n(l/2)n(l/2) so that n(l) oc exp{c/} up to power corrections. The constant c will depend on 
the detailed physics of the flux tube. Irrespective of its value, we see from eqn([TJ that there 
is a temperature 

T = T C = — (2) 

c 

at which the entropy compensates the Boltzman suppression, so that strings of arbitrary length 
are produced thermally. Clearly at T = T c the effective string tension, a e ff(T) = a — cT, 
vanishes and the system ceases to be linearly confining. 

This simple argument is illuminating in that it makes no reference to any underlying mi- 
croscopic theory, e.g. to gluons, and shows that any (linearly) confining theory must deconfine 
at some finite T. Since a e ff(T) — > as T — > T c , the corresponding correlation length diverges 
and we expect this phase transition to be second order. 

Given a phase transition at T = T c the partition function 

Z=J2 e ~ § ( 3 ) 

states 

should have some singularity there. Since the states contributing to Z are composed of 
glueballs, it is interesting to ask if there is any argument that Z should be singular at precisely 
the temperature T = T c = a/c at which confining strings 'condense into the vacuum'. We 
give the following simple argument. A natural model for glueballs is that they are composed 
of closed loops of confining flux. (See [S] for explicit models of this kind.) For a loop of length 
I the mass would be M ~ al. While this dynamics need not be convincing for small / and 
M (where the size of the state would be comparable to the width of the flux tube) it is hard 
to imagine excluding its contribution to the very massive glueballs that correspond to large /. 
Now for large / the number of states should be proportional to the number of different closed 
loops of length I. This number is clearly oc exp{c/}, with the same exponent c as in eqn(|T|). 
Thus Z will be singular at precisely T = T c = a/c with a Hagedorn-style divergence. (This 
classical number of states counting should be reliable for the relevant very heavy glueballs 
on the basis of the standard classical-quantum correspondence for large quantum numbers.) 
Note that the fact that the lightest glueball is very much heavier than T c , as are the splittings 
between the lightest glueballs, should not be regarded as paradoxical (see e.g. [§]); what 
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matters is the exponential growth of the number of glueball 'resonances' and the scale of this 
growth is, as we have just seen, precisely the same as in the earlier argument for deconfinement 
through string condensation. 

The above discussion has been formulated entirely within the effective low-energy confining 
theory. This means it does not depend on details, and so possesses an attractive universality, 
but it also means that it is necessarily incomplete. In particular it ignores the fact that in 
addition to the energy eigenstates composed of colour singlet glueballs, there are also states 
composed of gluons where the gluons are densely packed throughout the volume so that they 
avoid the constraints of confinement (in the same way that dense quark matter avoids these 
constraints). Neglecting any ambiguity in the separation between these two classes of states, 
we can write the partition function as 

Z ~ > e T + /_, eT= e T-\- e T (4) 

glueballs gluons 

where the Fq, F g are the corresponding free energies. Gluons clearly have an entropy S g oc 
0(N 2 ) in contrast to the Sg oc 0(N°) entropy of the colourless glueball states. (This is not 
quite correct because highly excited glueballs can be composed not only of closed loops of 
fundamental flux, but also out of closed fc-strings with k < N/2. But this is not important to 
the present argument.) This means that F g = E g —TS g will become equal to F G = E G — TSq 
at some T = T c and beyond T c the gluon states will be overwhelmingly more probable - we 
will be in the 'gluon plasma' phase. If this happens at a lower value of T than the 'string 
condensation' temperature discussed earlier, the deconfining phase transition will in fact go 
directly from a phase with a finite effective string tension to the non-confining 'gluon plasma' 
- and so we would expect it to be first order. 

The reason that the confining phase exists at all is that the confining vacuum energy is 
below that of the 'perturbative' vacuum (which one can roughly think of as the vacuum for 
the states that form the gluon plasma). This difference is often parameterised in terms of the 
gluon condensate and is usually expected to be ~ — 0(N 2 ). Thus low-lying hadron states are 
lighter by 0(N 2 )V than low- lying gluon states. Since at low T energy counts for everything 
and entropy for nothing, the gluon states are clearly irrelevant there. At large N the value of 
Fq will be dominated by this vacuum energy at all T, since all other pieces are 0(N°). So 
it is the competition between this and the —0(N 2 )T contribution to the gluon free energy 
F g that will determine T c . Since the vacuum energy density is on the scale of —0(N 2 a 2 ) we 
would expect 

lim T c = 0(» (5) 

N^oo 

and the fact that this is what one observes P] |2| may thus be regarded as evidence that the 
confining vacuum energy density is indeed —0(N 2 ). 

We know that for SU(iV > 3) gauge theories the phase transition is first order [TJ |2] and 
becomes more strongly so as N increases. This tells us that for these cases deconfinement 
occurs before the second order string condensation transition has had an opportunity to occur. 
For SU(2), however, the transition is indeed second order and the effective string tension 
vanishes as T — ► T~. It is therefore possible that here the phase transition proceeds by string 
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condensation. The same is true of both SU(2) and SU(3) in D = 2 + 1. In all these cases the 
deconfined phase immediately above T c might exhibit properties very different from those of 
a gluon plasma and more characteristic of what one would expect in the deconfined phase of 
a fundamental string theory - an interesting possibility whose investigation lies beyond the 
scope of the present paper. We note that if we assume the simplest effective string action, 
the Nambu-Gotto action, and if we assume it to be valid at all scales, then the deconfinement 
temperature occurs at [7j 



T c / 3 f 0.691 D = 3 + l , R . 

(6) 



yja y 7Tp-2) t 0.977 D = 2 + 1 
which we can compare to the lattice values [21 El El : 



rp ( 0.709(4) SU(2), D = 3 + 1 

1.121(8) SU(2), D = 2 + 1 (7) 
0.985(12) SU(3), D = 2 + l. 



These predicted and lattice values are remarkably close (as has been previously observed 
jHl), simultaneously suggesting that the Nambu-Gotto action provides a surprisingly accurate 
representation of the dynamics of the confining flux tube all the way down to scales of 0(T C ) 
and that these phase transitions are indeed driven by 'string condensation'. 



2.2 Multiple transitions? 

We have assumed in the above discussion that there is a single deconfinement transition. 
However this need not be the case. As is well known, deconfinement is associated with the 
complete spontaneous breaking of a Zn symmetry, and for N > 4 this symmetry could be 
broken in more than one step. In the case of SU(4), for example, one could imagine a breaking 
Z 4 — > Z 2 at T = T c followed by a breaking of the remaining Z 2 at T = T c > > T c . Such a 
pattern of breakings occurs in certain spin models ^Uj and so is not entirely idle speculation. 

Physically such a step-by-step breaking of the Z N symmetry can be related to the decon- 
finement of different /c-strings at different temperatures T c (k). Recall (see e.g. [II]) that at 
low T a source that is like a local clump of k fundamental charges will be joined to a distant 
conjugate source by a stable flux tube that is called a fc-string (for k < N/2). One can think 
of a fc-string as a bound state of k fundamental strings. The tension of such a string is found 
to lie somewhere between the Casimir scaling 12J and 'MQCD' ^3] conjectures [TH I14j: 

- fc) < £fc < s j^w_ (8) 
N- 1 - a ~ sin f 1 ' 

The two-step deconfinement that we outlined in the previous paragraph would correspond to 
k = 2 charges becoming deconfined at T = T c while k — 1 charges remain confined until a 
phase transition at a higher temperature T = T c i. This can be generalised in an obvious way 
to larger iV and k. 
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When deconfinement proceeds via string condensation, it is easy to see what is needed 
for the above scenario to occur. The entropy factor in eqn(£Q) will in general depend on k, 
nkil) oc expcfcZ, so that fc-strings will condense at 



If we had no knowledge of the detailed dynamics we could easily imagine that the Ck might 
be such that a T c (k = 2) < T c (k = 1). (Note that the reverse is not possible because once 
k — 1 sources deconfine, all k > 1 sources will also deconfine.) However the effective string 
theory arises from a gauge theory and we will now argue that in a gauge theory this probably 
cannot occur. 

First we note that Ck must contain a factor 1/4 that provides a scale for the / in the 
exponent. This scale If. will be given, in the simplest case, by the radius, r^, of the fc-string. 
This assumes that different flux tubes are geometrically much the same apart from their radii. 
Now, a flux tube carries non-Abelian electric flux and if we assume that the flux is distributed 
uniformly across the cross-section of the flux tube, a simple classical calculation would predict 



where we have used the fact that the trace is the quadratic Casimir and we have inserted the 
value appropriate to a fc-source in the totally antisymmetric representation. (This classical 
calculation is just like the old bag model ^H] but without the variational determination of the 
radius.) Inserting this and Ck = co/r^ into eqnflnj) we find 



Although the argument is clearly too simple, the ratio is so much larger than unity that 
it seems very plausible that fundamental sources must deconfine at the lowest temperature, 
leading to a simultaneous deconfinement of all other sources. 

For N > 3 and for fundamental sources the phase transition is first order, passing directly 
to a 'gluon plasma' without string condensation. This dynamics will presumably break the 
centre symmetry completely so that there is a single T c . The possibility that there might 
be some earlier second order phase transition corresponding to the condensation of higher k- 
strings is unlikely given our arguments above coupled to the fact that the first order transition 
temperature is less than the string condensation temperature of fundamental strings. 

Although our conclusion is that there is almost certainly a single deconfinement transition, 
the possibility that things are otherwise is too interesting to ignore. We therefore perform a 
numerical investigation in Section although without unearthing any surprises. 

2.3 V and N dependence 

When the deconfining phase transition is first order, one determines T c as lying in the range 
of T where the system undergoes back-and-forth tunnelling between the confined and decon- 



T c (k) 



(9) 




(10) 
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fined phases. This tunnelling necessarily proceeds through intermediate field configurations 
containing two interfaces that divide the periodic spatial volume into roughly two equal parts. 
Thus, at T = T c where the two phases have the same free energies, the tunnelling probability 
will be suppressed as 

Ptunnel OC exp{ -2(T cd A/T c } (12) 

where a c d is the surface tension of the interface and A is the (smallest) spatial cross-section, e.g. 

2 

A = V$ for a symmetric volume V . So to be able to locate the actual phase transition we need 
to work on a small enough volume for this tunnelling to be non- negligible. In that case, such 
back-and-forth tunnelling will occur for a range AT of temperatures around T c for which the 
difference in the free energies AF/T is 0(1). (Once AF/T ^> 1 the transition will effectively 
just go one way.) For T close to T c we can use the approximation AF ~ L h V(l — T/T c ) where 
Lh is the latent heat (per unit volume) of the transition, and we see that AT and hence the 
uncertainty in locating T c on the finite volume, is given by 

AT cc -2- (13) 

Thus we see from eqn fllHj) that the latent heat determines the uncertainty in determining 
T c on a volume V, while eqn ()12)) tells us that the interface tension determines how large a V 
one can consider and still see the transition at T = T c . Of course, even if V is so large that no 
back-and-forth tunnelling can be seen at T = T c , if we move T away from T c then the growing 
free energy difference will eventually drive the tunnelling in one direction. The tunnelling will 
proceed by the formation of a small bubble that grows, improbably, to a certain critical size 
beyond which the free energy difference between the two phases outweighs the cost of the 
bubble interface free energy. At this point the bubble will rapidly grow to fill the available 
volume. (Such a continuous description is appropriate since the lattice Monte Carlo proceeds 
by changing one link matrix at a time, and so is effectively a continuous deformation of the 
fields.) It is clear that as the bubble grows, the cost of the bubble will first grow and then 
decrease. If we estimate the tunnelling probability by the thermodynamic suppression of the 
most suppressed bubble in this process, we find that for 

(T - T c ) 2 > 4fT c (14) 

(c is some numerical constant) one-way tunnelling will occur however large is the volume. 
This gives the width of the hysteresis and provides a finite uncertainty on our estimate of T c if 
we work on very large volumes. To calculate T c one therefore works on volumes that are small 
enough to possess tunnelling right across the transition and one then extrapolates to infinite 
volume. On the other hand, if we choose to work with V so large that no tunnelling occurs in 
a (realistic) simulation at T ~ T c , we will be able to calculate properties of the confining (or 
deconfined) phase at T = T c - and in its neighbourhood. We shall make use of this strategy 
for N > A where the transition is strongly enough first order for this to be practical. 
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In an earlier paper [2] we performed a calculation of the latent heat at a fixed value of the 
lattice spacing, a ~ 1/5T C , and found that 

L h N oT N\ (15) 

just as one expects from the usual large- N counting arguments. In the present paper we will 
calculate the continuum limit of for N = 4, 6, 8 and so determine not only the variation 
with N but also its value in physical units. In [2] we also showed that the interface tension 
grows with N. In the present paper we attempt to distinguish between o c d oc N and a C( i oc N 2 
using a different method, but unfortunately our accuracy will not allow us to answer this 
question convincingly. (There is a simple theoretical motivation for a c d oc iV 2 that has been 
given in [2].) Each of these two possibilities leads to a particularly simple picture of the phase 
transition as iV — > oo. We first note from eqn(|12Jl that for either power of N the back-and- 
forth tunnelling when T = T c and AF = will be exponentially suppressed in iV (or in N 2 ) 
however small we make V. However we also see from eqn()15|) that if a c d oc N then the width 
of the hysteresis will be no more than AT oc 1/N. Thus as N — > oo on a fixed volume V, 
however small or large, the phase transition will become completely sharp. If, on the other 
hand, cr cd oc iV 2 then this is no longer the case. In this case egn lfTlj) tells us that on a large 
volume the hysteresis grows as AT oc iV i.e. becomes infinite at iV = oo. In this limit of 
infinite hysteresis the confined and deconfined phases exist at all T with no tunnelling between 
them. (If V is smaller, below the critical bubble size, the contribution of the interface tension 
always dominates as the bubble grows, so this infinite hysteresis scenario would in fact occur 
on all volumes.) Of course our argument for eqn ()14j) breaks down once AT becomes large, and 
in any case the string condensation transition will, presumably, intervene somewhere above 
T 

These two ideal scenarios - zero or 'infinite' hysteresis on any volume V - provide one 
motivation for being interested in the value and iV-dependence of the interface tension a c d- 
Unfortunately the numerical results we present in Section |H] will not be conclusive. All this 
is also reason to be interested in the latent heat although, since that is a fundamental quan- 
tity characterising the phase transition, we would be interested in it in any case. Here our 
numerical results, in ,2] and Section are more precise. 

2.4 Masses 

The size of the latent heat (rescaled by N 2 ) tells us how strongly first order the phase transition 
is. Another measure of this strength is the mass gap in each of the phases as T — > T c . For a 
second order transition the mass gap vanishes; for a weakly first order transition it becomes 
small but non-zero; for a normal first order transition it should be on the order of the typical 
dynamical length scale of the theory, e.g. T c or the T = string tension. If the mass gap 
increases with N it provides us with additional evidence that the transition is becoming more 
strongly first order. We provide some calculations of this mass gap in Section |3 

Since the transition is to do with deconfinement, the natural mass to look at is the effective 
string tension at finite T, a e ff(T). As discussed, for example, in Section 3.1 of j2] this is given 
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by the mass m t (lt) of the lightest flux loop that winds around the timelike torus: 



°eff{T) = jm t (l t ) = Tm t {l t = i). (16) 

Thus at T = T c it will be this mass, m t (l t = l/T c ), that will vanish for a second order 
transition, and should be ~ T c for a normal first order transition. 

Note also that egnffTB^) allows us to derive the dependence of a e ff(T) on T for low tem- 
peratures, T ~ 0, in a very simple way. We recall that the leading l t correction to the linear 
dependence of m t {l t = 1/T) is provided by the universal bosonic string correction, n/3l t , so 
that eqn (fi7)|) implies 

a eff (T) T =°a-^T 2 + 0(T% (17) 

A parallel argument holds for /c-strings and their associated string tensions, at any fixed N. 
If we take N — > oo at fixed k, the /c-string becomes an ever more loosely bound state of k 
fundamental strings and, for reasons discussed in Section 3.1.2 of [Hj, this implies that the 
leading correction in eqn ()17j) will hold in an ever decreasing interval around T = 0. 

This winding flux loop is the lightest state coupling to timelike Wilson lines (Polyakov 
loops). In the deconfined phase this operator acquires a non-zero vacuum expectation value. 
At high T it can be expanded and the first non-trivial term is oc Tr{A$} showing that 
the strongest coupling of the vacuum subtracted operator will be to two (electric) Debye- 
screened gluons. To leading order in the 't Hooft coupling, the Debye mass must clearly be 
m 2 D oc g 2 (T)NT 2 on dimensional grounds, and a one- loop calculation ^Hj shows that in fact 



m t (T) = 2m D = 2yJ 9 ^ N T. (18) 

Near T = T c this perturbative calculation and interpretation is no longer compelling but we 
shall use the same language for simplicity. 

A quite different quantity from the flux loop that winds around the timelike torus is the 
one that winds around a spatial torus. As T — > this difference disappears of course and 
either will provide a measure of the T = string tension. At high T the mass, m s (l s ), of this 
spacelike flux loop and its associated string tension, a s , will vary as 

cr s (l s ;T) = T m s (l s ) oc T 2 (19) 

on dimensional grounds. The 'spacelike' string tension, a s , becomes the string tension of 
the linearly confining D=2+l gauge-scalar theory to which the D=3+l SU(iV) gauge theory 
is dimensionally reduced at high enough T. In principle the spatial string tensions in the 
confined and deconfined phases at T ~ T c could be the same, and this has often been inferred 
from earlier calculations in SU(3). However when one works at higher iV where the phase 
transition is quite strongly first order, so that the separation into the two phases can be easily 
made unambiguous, one can show, as we shall in Section |U1 that there is in fact a significant 
discontinuity in a s across T c . 
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2.5 Thermal average or Euclidean artifact? 



In using the Euclidean time framework for calculating thermal averages at a temperature T 
it is important to be aware that dynamical features of the 4 dimensional Euclidean fields will 
only be features of the finite temperature field theory if they are encoded as such in the (equal 
time) thermal averages. 

A cautionary example concerns the well-known spontaneous symmetry breaking of 
SU(iV) gauge theories for T > T c . It is certainly the case that as the Euclidean time extent, 
It, of our 4-volume decreases through the critical value It = 1/T C , the 4-dimensional system 
undergoes a phase transition to one of N degenerate vacua. The timelike Wilson line acquires 
a value that is proportional to one of the iV'th roots of unity and this serves as an order 
parameter distinguishing these various phases. There are interfaces between these vacua and 
these are perturbatively calculable at high enough T [T7J [TH]. Indeed we note that if we 
relabel co-ordinates (x, y, z, t) — > (t, x, y, z) so that the short torus is now in the ^-direction 
and the Euclidean-time torus is now long, so that the theory is effectively at zero temperature, 
then this phase transition becomes a real finite volume phase transition, obtained by reducing 
l z while keeping the other three lengths fixed (and large). All the properties listed above 
(suitably relabelled) become physically observable. Nonetheless, despite appearances, it is 
now widely accepted that all this does not mean that the deconfined phase of SU(iV) gauge 
theories possesses such a symmetry breaking with, for example, bubbles of different Zn vacua 
in the high-T gluon plasma. At least, no-one has succeeded in showing that such a structure 
is encoded in an appropriate fashion in equal-time thermal averages. (The timelike Polyakov 
loop clearly cannot serve such a purpose.) 

It is therefore always important to consider if what one is seeing might be a Euclidean 
artifact rather than a genuine property of the gauge theory at finite T. This applies to the 
details of the phase transition itself. Of course the first order transitions for SU(iV > 3) are 
reflected in discontinuities in thermal averages and there is no ambiguity. However the second 
order SU(2) transition is more problematic. In particular the diverging correlation length - 
easily seen in the correlation of timelike Wilson lines - is in practice not visible in equal-time 
correlators and one finds no growing peak in the specific heat. It is therefore interesting to 
ask in what way this transition manifests itself in the high-T gauge theory, and this we do in 
Section El 

3 The lattice setup 

Our Euclidean space-time is discretised to a hypercubic, periodic lattice of size L\L t and 
lattice spacing a, with lengths denoted by L when there is no ambiguity. We assign SU(iV) 
matrices, Ui, to the links I (also denoted by U^n) for the link emanating in the [i direction 
from site n). We use the standard plaquette action 




(20) 



v 
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where U p is the ordered product of the SU(iV) matrices around the boundary of the plaquette 
p. (We shall often use the shorthand u p = ReTrU p /N.) S appears in the Euclidean Path 
Integral as 

Z= ! Y[dU ie - s (21) 
J i 

and becomes, in the continuum limit, the usual Yang-Mills action with 

2N 

P=— (22) 
9 

Since the theory is asymptotically free, we approach the continuum limit by reducing the bare 
coupling, g 2 (a) — ^ 0, and so (3 —>■ oo. Our simulations in this paper are performed with a 
combination of heat bath and over-relaxation updates in which we update all the N(N — l)/2 
SU(2) subgroups of the SU(iV) link matrices, as described in (THJIH]. 

The thermodynamic partition function of the gauge theory at a temperature T is identical 
to the Euclidean Feynman Path Integral with temporal periodicity of 1/T. So by calculating 
appropriate expectation values on an L\L t lattice with L s 3> L t we are calculating a lattice 
approximation to the corresponding thermal average at 

(23) 



which becomes exact when we extrapolate to the continuum limit at fixed T. 

When we vary iV we expect |2lH l2~Tj that we will need to keep constant the 't Hooft coupling, 

A, 

A = 9 2 N (24) 

to obtain a smooth large- N limit. Non-perturbative calculations for 2 < iV < 5 have supported 
this expectation Using the string tension calculations in our recent paper ^3] we can 

extend this study to 2 < N < 8. Following [TH] we define a (lattice improved) 't Hooft 
coupling on the scale a by 

2N 2 

Xj(a) = 9l (a) 2 N = —— (25) 

P{Up) 

and plot in Fig^the result. We observe very nice evidence that the running 't Hooft coupling 
approaches a limiting iV-independent function of the length scale I = a^/a as iV — > oo, and 
that it does so very rapidly, once we are away from the coarsest values of a that are close to 
the bulk transition/cross-over (see Section ITUJl . 
Finally two comments for the general reader. 

Our Monte Carlo generates a sequence of fields by changing the value of one link matrix at 
a time, with the local action density controlling the size of the change. Thus the Monte Carlo 
proceeds by what is essentially a local deformation of the fields and so it mimics features of a 
real phase transition in a statistical mechanics system. That is to say, a first order transition 
can be described as a continuous tunnelling process, in terms of growing bubbles of one phase 
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within another phase, and if at T = T c that probability is low enough then we will observe a 
hysteresis effect. 

Our second comment is that the deconfining phase transition occurs only on physical length 
scales. Thus in lattice units the latent heat will be a 4 Lh oc a 4 and the mass gap at the 
transition will also satisfy am g 0. So from the point of view of the lattice, the transition 
is very weakly first order and becomes second order as a — > 0. This is however an illusion, 
merely reflecting the fact that the continuum limit is a second order transition of the lattice 
theory. When we express the mass gap and latent heat in physical units, e.g. as L^/T 4 and 
m g /T c , it becomes clear that the transition is robustly first order for N > A. 



4 The deconfining temperature 

Since presenting our results on T c /y/a in [2j we have significantly improved our SU(8) calcu- 
lations. We now have a greater range of volumes at L t = 5 and this allows us to obtain a 
more reliable estimate of the coefficient h of the leading 1/V correction to f3 c : 

P c (V) V =°°p c (oo)-h^. (26) 

We also have better calculations of the SU(8) lattice string tensions [Hj and this allows us to 
perform a more accurate continuum extrapolation of T c /^Ja. 

We recall [2] that our strategy is to use the value of h calculated at a = 1/5T C to perform 

V — > oo extrapolations at the smaller values of a at which we usually have performed cal- 
culations on only one volume. In j2] we did not attempt to compute the a-dependence of h 
and we simply doubled the statistical error on h in the hope that this would encompass any 
variation of h with a(/3). However the form of this a dependence of h is known |2| 

h ( a ) a =° h + (a 2 ) (27) 

and so we can use the values of d\na/d(3 that we calculate later on in this paper to rescale h 
at the P values of interest. We shall do this not just for N = 8, but for N = 4, 6 as well. We 
then increase the error on h by 15% in order to cover the remaining 0(a 2 ) corrections which, 
given what we observe with other dimensionless quantities, we expect to be small. 

In Tables [U - El we list our calculated values of the critical /3, obtained by extrapolating to 

V = oo using eqns ()26l27|) with the values of h as shown. We calculate the string tension at 
each f3 using the values in and this provides us with the listed values of T c j y/a. Comparing 
with Table 5 of [2] we see that the SU(4) and SU(6) values have changed very little, while the 
SU(8) values, where we employ better calculations of a, have changed more substantially. 

We extrapolate to the continuum limit using an a 2 a correction, as described in and list 
the resulting continuum values of T c /^/a in Table |U For completeness we have included the 
SU(2) and SU(3) values although their analysis is unchanged from [2]. We note that the fits 
all have an acceptable x 2 P er degree of freedom except for SU(4) where the fit is very poor. 
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This promises to create a problem for any statistical analysis that involves the SU(4) value, 
as for example in the extrapolation to N = oo to which we now turn. 

In FigEJwe plot our continuum values of T c j \/a against 1/N 2 . The usual large- N counting 
tells us that the leading large- N correction should be 0(1/N 2 ), and this corresponds to a 
straight line on this plot. While it is clear that our results do fall approximately onto a 
straight line, the best fit \ 2 ls m f ac t very poor. It is reasonable to suppose that this is due to 
our SU(4) value being off perhaps by several standard deviations. To try and deal with this 
we rescale the errors on the SU(4) lattice calculations by almost a factor of two, so that the 
continuum extrapolation has a x 2 P er degree of freedom of unity. This increases the error on 
the continuum extrapolation to the value shown in square brackets in Table EJ Performing a 
fit with this modified SU(4) error, we obtain 

^ = 0.5970(38) + °- 4 t 9 2 (29) (28) 

with an acceptable \ 2 P er degree of freedom of about unity. This fit is plotted in Fig|21 As 
we have noted before PjJ |2] it is remarkable and perhaps puzzling that the leading correction 
suffices all the way down to SU(2) despite the fact that the transition changes from first to 
second order between N = 3 and N = 2. 

In performing our large- N extrapolation we assume a conventional 0(1/N 2 ) correction. 
Since this expectation is based on diagrams (albeit to all orders) while T c j^Jo is a non- 
perturbative quantity, it is interesting to see whether our results are accurate enough to 
constrain the power of the leading correction. We therefore perform fits to a 1/N a correction. 
We find that the \ 2 P er degree of freedom of the best fit is 



3.5 
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= 1 


0.8 


a 


= 2 


3.1 
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= 3 



(29) 



*!.= 

which provides some evidence that the leading correction is indeed 0(1/N 2 ). 

5 Latent heat 
5.1 Background 

The free energy, F, and the free energy per unit volume, /, are defined to be 

EE F fV , 

e ~T = e ~T = e -— (30) 

states 

where V is the spatial volume and Z is the standard finite temperature partition function. 
The internal energy density, e, is the average energy per unit volume 
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Suppose we are on a L\L t lattice that remains fixed in lattice units. The Euclidean 
Feynman Path integral provides a definition of a thermodynamic partition function of the 
lattice system at T = l/aL t which will tend to the desired continuum partition function as 
we reduce a (as long as the spatial volume is sufficiently large), a and hence T are varied by 
varying (3, so eqn lpTTj) becomes 

T 2 d(3d T 2 r 3 ^ (tip) d(3 

where (u p ) is the average plaquette. We thus see that the latent heat in units of T c is given 
by 

^ = ^ = -L t a-6A(u p ). (33) 

Note that if we wished to calculate the internal energy or the free energy, we would need to 
be concerned about regularising it, for example by subtracting from (u p ) in eqn (}32|) the value 
of the plaquette at some reference value of T (but at the same value of a) . Since we are only 
interested in the difference in internal energies, such issues need not concern us here. 

The finite T phase transition is clearly also a phase transition of a statistical system in 4 
space dimensions whose classical partition function is our Euclidean Feynman path integral. 
In this reinterpret at ion, f3 is the inverse temperature and the action S is the energy. In this 
case the latent heat is simply the jump in the action i.e. it is just what we find in eqn (|33j) 
but without the factor of ad (3 / da. This latter factor is, however, crucial if we want to obtain 
an estimate of the real deconfining latent heat of the gauge theory in 3 space dimensions. 

We can calculate adf3/da using 

a d ± = ail J^- (34) 
da d(afi) 

where \i is some physical quantity with dimensions of mass, e.g. /i = ^Ja or fi = T c . Such 
different choices will differ by 0(a 2 ) lattice corrections in Ae/T c 4 . We shall choose to use 
fi = ^Jo . To obtain the value of d(a^a)/d(3 at the values of j3 that correspond to T = T c for 
various values of L t we interpolate between our calculated values of a^/a (as listed in Tables 
5-9 of |I3j) using a fitting function 

W) = Coexp{-^ i ^j^ + c 1 x + c 2 x 2 + c 3 x 3 } ; x = iV 2 (l - 1). (35) 

The value of (3q is chosen to coincide with one of the calculations near the middle of the 
interesting range of (3. (Given that we fit cq, this is not an extra parameter.) Where we can 
obtain a good fit without the 0(x 3 ) term we do so. Note that although the fitting function is 
designed so that it automatically tends to one-loop scaling at large enough /3, the particular 
form of the fit does not really matter as long as we intend to use it to interpolate rather than 
to extrapolate in (3. (Including a In (3 term corresponding to 2-loop scaling does not improve 
the fits and does not even appear to improve the rate at which scaling is approached when the 
fits are extrapolated to larger (3.) In Table we list the values of the parameters of the best 
fit for each N together with the range of j3 interpolated (and the y 2 per degree of freedom of 
the best fit). 
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5.2 Results 



In j2] we calculated the value of A(u p ) from the V — > oo limit of the specific heat peak. 
(Note that in [2j we referred to this as Lh, and the latent heat was defined as the jump in 
the average 'energy' per plaquette of the four dimensional lattice system. This differs from 
the physical latent heat of the three dimensional thermal gauge theory, as discussed above.) 
This calculation requires a detailed finite volume study, which we have only carried out at 
a ~ 1/5T C , with the results shown in Table |HJ and so it is only at this fixed lattice spacing 
that we have been able to compare latent heats for different N. What we would now like to 
obtain is the comparison in the continuum limit, and for that we need a different technique. 

We begin by observing that on our very largest lattice volumes at a ~ 1/5T C there is no 
tunnelling for iV > 4 at T ~ T c in our simulations. We can therefore produce ensembles of 
lattice fields that are entirely confining or deconfining by obtaining the starting lattice field 
from T — > T~ and T — > T c + respectively (and subsequently thermalising at T = T c ). From 
these two ensembles we can calculate (u p ) c and (u p )d, respectively, and hence A{u p ). The 
values of the latter are listed in the L t = 5 rows of Table [H They are all compatible with 
the values obtained from the specific heat extrapolation, as listed in Table reinforcing our 
confidence in both methods. 

We can extend this method to smaller volumes where the tunnelling is present but is 
relatively rare. We can use the mean Polyakov loop as an order parameter to categorise sub- 
sequences in the Monte Carlo generated sequence of lattice fields as being confined, deconfined 
or (much more rarely) tunnelling. There is clearly some ambiguity in this separation, but this 
ambiguity disappears as the tunnelling becomes less frequent. We try to incorporate this am- 
biguity into the error estimate by using various different cuts on the order parameter and by 
averaging the order parameter over subsequences of lattice fields of different length. However 
once the spatial volume becomes sufficiently small in physical units - as measured by L s /L t 
- all these remedies become inadequate and our values will be increasingly unreliable. A test 
of this method for the case of SU(4) can be found in Fig. 5 of pQ where the L = 32 lattice 
has no tunnelling and provides the benchmark for the comparison. We observe that the error 
becomes rapidly larger as we decrease V, as it should, but that the value of A{u p ) remains 
consistent with the V — > oo value, even on the smallest volumes where some separation be- 
tween phases is apparent. This reassures us that this method can be used at smaller values 
of a where we do not have calculations on volumes that are very large in physical units. 

Our results are listed in Table [7| For N > 4 we have three values of L t . Since our 
calculations are at a fixed temperature, T ~ T c , we have a oc 1/L t and this will allow us to 
perform a continuum extrapolation. At each N and L t we have used only the largest lattices. 
We show the value of (3 at which A{u p ) is calculated. In fact around each of these values of 
(3 there are usually two or three nearby values which exhibit tunnelling and are useful. We 
transport these values to the indicated value, varying a 4 using the listed values of 8(3/ din. a. 
We then transport from these reference values of (3 to the values of /3 C (V = oo; L t ) listed in 
Tables [OH] and using eqn (J33|) we obtain the latent heat at T = T c in physical units, as shown. 
We now discuss these results. 

Our best values of Lh, in the sense of possessing the smallest statistical and systematic 
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errors, are those for L t = 5. In FigJ3]we plot these latent heats against 1/N 2 and show a fit 
of the form 

^{^} 4=Co + ^- (36) 

This form is expected because the average energy of the gluon plasma should be 0(N 2 ) and 

the leading large- N correction should be 0{1/N 2 ). We see that L h /N 2 does indeed appear to 

go to a finite limit and that the SU(3) value is well below the fit, in line with the usual belief 

that it is weakly first order. We note that Lh/N 2 increases to an N = oo limit that is 0(1) in 

units of T c . That is to say, the transition is robustly first-order in that limit. 

i 

In Fig|3]we plot the dimensionless ratio L^/Tc against a 2 T 2 . For the plaquette action the 
leading lattice correction is known to be 0(a 2 ) so we expect to be able to extrapolate to the 
continuum limit using 



s + ca 2 T 2 (a). (37) 
T c (a) T c (a = 0) cK J K J 

Unfortunately it is clear that such straight-line fits will be very poor and indeed they lead to 
X 2 / n df — 3,9, 11 for N = 4, 6,8 respectively, which is unacceptable. All the plots show the 
similar feature that the L t = 5 value is much too high to lie on the fit. Although in some cases 
L g / L t is smaller for L t = 6 than for L t = 8, in other cases it is larger, so we do not attribute the 
problem to large finite volume corrections afflicting all the L t = 6 calculations. We therefore 
take this to be a real effect rather than a manifestation of some unexpectedly large systematic 
error. That means that we must extrapolate using just the two values, at L t = 6 and at 
L t = 8. This means we have no evidence that including just the 0(a 2 ) correction is sufficient. 
Nonetheless, it is plausible to assume that it is since for other quantities, such as glueball 
masses, the deconfining temperature, etc., we know that it is in this range of a. Making this 
assumption we obtain the continuum values listed in Table |Hl using the extrapolations shown 
in FigH 

i 

We can now take these continuum values of L^/T c and see whether they determine the 
N = oo value. We find that they are indeed well fit by the expected functional form 

4 - 0.766(40)-!^ : = 0* (38) 



NiT c ' v ' N 2 

although the errors are rather large because of the fact that we could not use our most accurate 
data, which was at L t = 5. 

As we shall discuss in Section [TUJ the value of (3 corresponding to a = 1/5T C is very close 
to the bulk transition. This is a rapid crossover for N = 4, and becomes a strong first order 
transition for N > 6. Indeed, for N = 8, a = 1/5T C is very close to the lower f3 edge of the 
bulk hysteresis curve. At this first order bulk transition the jump in the plaquette, A{u p ), is 
about 50 times larger than it is for the deconfining transition on the L t = 5 lattice system. 
It is therefore plausible that the whole effective potential is becoming deformed at this value 
of a and that this is causing the anomalously large value for the latent heat at L t = 5. The 
fact that the effect is larger for N = 6,8 than for N = 4, where the transition is only a 
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crossover, supports this hypothesis. This transition can be naturally regarded as a finite N 
manifestation of the D=l+1 Gross- Witten transition [221 12H] and the plaquette is the natural 
order parameter. It is therefore perhaps plausible that it should be the average plaquette that 
is most sensitive to the nearby presence of the bulk phase. 

6 Masses 

At T = T c we can perform calculations separately in the confined and deconfined phases by 
working with volumes large enough for there to be no sign of any tunnelling (or attempted 
tunnelling). We have performed such calculations on 32 3 5, 16 3 5 and 12 3 5 lattices for SU(4), 
SU(6) and SU(8) respectively. In SU(3) we use a 64 3 5 lattice, on which there is tunnelling, 
but it is very rare and we can work with long subsequences of field configurations that are in 
a single phase. (The reason that we are able to work on smaller spatial volumes, V = Z 3 , as 
N increases is because in the tunnelling probability, oc exp{— 2<r cd l 2 /T}, the interface tension, 
a c d, increases with N - as discussed in Section EJ) Although the calculations in this Section 
are performed at a fixed value of the lattice spacing, this value is small enough that the results 
should be representative of what one would obtain in the continuum limit. 

6.1 Mass gap at T = T c 

We calculate the lightest mass that couples to the Polyakov loop which winds around the 
time-torus. We do this for different fc-loops (see Section 12. 2|) at each value of N and we 
present the masses for the confined and deconfined phases in Tables 19*1 and I1UI respectively (In 
the deconfined phase we use the vacuum subtracted loop.) 

We observe that in the confined phase the higher k loops are more massive, just as they are 
at T = 0. If we take the mass of the SU(8) loop in the confined phase, and apply the leading 
string correction (a priori inadequate for such a short string) so as to extract a string tension 
a = m/l t + n/3l 2 , we obtain, for SU(8), a value, ay/a ~ 0.314(3) that is quite close to the 
T = value ay/a ~ 0.346 [Tff|. This suggests that this first order transition is quite strong. 
If we express the loop mass in the natural units of T c and plot it, in Fig|"51 against 1/N 2 we 
see that this mass gap tends to a finite N = oo limit, albeit with a large 0(1/N 2 ) correction. 
Indeed we see that for SU(3) the correlation length, £ t = l/m t , is quite large indicating a 
transition that is quite weakly first order. 

The k = 2 mass gap in Fig "5] also appears to tend to a finite N = oo limit. Although 
we have used a leading 0{1/N 2 ) correction, this is not guaranteed to be correct for k > 1 
strings. The alternative is an 0(1/ N) correction which arises in, for example, the Casimir 
Scaling conjecture f""*"| . It is interesting to note that the string tension ratio 

o- ef f(k = 2) m t (k = 2) f 1.91 ± 0.13 5a k=2 = 0(1/N 2 ) 

a eff (k = 1) m t (k = 1) ~~ * 1 2.36 ± 0.21 5a k=2 = 0(1/ N) [6J) 

is compatible with the value of 2 that one expects from factorisation in the N = oo limit (even 
at finite T) for both forms of the correction. 
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In the deconfined phase all the k- loops couple to the same object (the centre symmetry is 
spontaneously broken) and, as we see in Table ITUl their masses are the same. (It is interesting 
to note that the operator with the largest k appears to have the best overlap.) Writing 
this mass as twice the Debye mass, we plot its value against 1/N 2 in Fig® Once again we 
see a finite large-iV limit with a large finite-iV correction that leads to a small SU(3) mass, 
confirming its weakly first order character. 

It is of course quite possible that some of our calculated masses suffer significant finite 
volume corrections. To partially address this question we show in Table ^2 the masses one 
gets in SU(4) on smaller lattices at the same value of a. The smaller the lattice the greater 
the tunnelling and the more difficult it is to distinguish the two phases, hence the increasing 
errors. To partially control this we also show some calculations performed at slightly different 
values of (3 where there is very little tunnelling. The conclusion of this comparison is that the 
confined k — 1 mass has at most modest corrections but that they may be more substantial for 
the k = 2 mass. In the deconfined phase, any finite V correction to the L s > 20 calculations 
look small. Since finite V corrections are likely to be functions of the lattice size in units 
of the mass, the above suggests that our SU(6) and SU(8) lattices are large enough for the 
confined k = 1 masses and probably large enough for the deconfined masses, but that the 
k > 1 confined masses may well suffer substantial corrections. (See also [24J.) 

6.2 Debye mass at T > T c 

We see from eqn(|18|) that at leading order in g 2 (T)N the Debye electric screening mass 
is independent of N and grows linearly in T up to the weakly varying coupling factor. 

We extract this mass from correlations of the vacuum-subtracted timelike Polyakov loop 
on L s = 8 lattices with L t = 2, 3, 4 at a value of (3 where L t = 5 corresponds to T = T c and list 
the results in Table IT21 Thus our lattice calculations cover the range 1.25 < T/T c < 2.5. In 
the case of SU(8) the tunnelling on an 8 3 5 lattice is sufficiently well-defined that we are able 
to separate fields in the confined and deconfined phase and so calculate mo down to T = T c . 
For T = T c we also have calculations on the larger volumes listed in Table [HI We see that 
at T = T c , where the masses are smallest, the finite V shift is at the level of ~ 20%. At the 
highest T, where the masses are largest, we have also performed calculations on 6 3 2 lattices, 
for SU(4) and SU(8), and there is no finite V shift visible within quite small errors. 

In FigJ7|we plot m t = 2mn against T/T c for both SU(8) and SU(3). We observe that for 
T > 5T c /3 the mass is independent of iV and appears to increase oc T. There are only small 
deviations from this behaviour at T = 5T c /4 and it is only at T ~ T c that we see any large 
deviation. 

Thus the asymptotic T — > oo behaviour sets in very close to T c , especially at larger N. 
This is remarkable since for T ~ T c the coupling will surely not be small. Indeed, if we apply 
eqnflTHJ) to our calculated values, we find 

9 2 (T^1.25T C )N ^ o 43 
An 

Such a substantial coupling also means that the lowest order calculation is probably unreliable 
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(as may be the interpretation of m t as 2m D ). Indeed eqn ffTSj) would predict that the approach 
to the linear behaviour at high T should be from above, since g 2 {T) decreases with increasing 
T, whereas what we observe in Fig0 is that it does so from below. 

6.3 Spatial string tensions at all T 

The expectation values of large spatial Wilson loops are known to decay exponentially with 
their area at all T and the coefficient of the area is the 'spatial string tension'. At high T this 
corresponds to the string tension of the dimensionally reduced D = 2 + 1 SU(iV) gauge-Higgs 
theory. It also corresponds to the string tension in a L 2 s L t spatial volume at a temperature 
aT = 1/L S , for sources separated along one of the long directions. (This follows from a trivial 
relabelling of the Euclidean co-ordinates.) 

The spatial string tension, <J S (T), can be obtained from the correlation of Polyakov loops 
that wind around a spatial torus. From this we obtain the 'mass', m s , of the lightest flux 
loop that winds around that torus and applying the usual string correction, which should be 
adequate since L s is large, we obtain the spatial string tension: 

^am s (L s ; T) = a 2 a s {T) - ^-L. (41) 

Here D± is the number of dimensions transverse to the Wilson loop. For D=3+l and low 
T we have D± = D — 2 = 2 and we would expect this to remain the case for any T in the 
confining phase. It is less clear what one should use in the deconfined phase: at very large T 
one presumably has = 1. Since eqn (|T9"]) tells us that cr s (T) oc T 2 , this correction will be 
0(1/T 2 ) and so unimportant at high T if L s is kept fixed. 

Below T c er s (T) is more-or-less constant and above T c it grows with T. Is the change smooth 
or discontinuous? Given that the transition is first-order a discontinuity would be natural but 
the existing folklore, based on SU(3) calculations, is that it is continuous. To address this 
question we show in Table IT3*1 the values of am s as obtained on sequences of 12 3 5 lattices at 
T = T c in SU(8), that are entirely in confined and deconfined phases respectively. (These 
volumes are large enough to suppress any sign of tunnelling.) What we actually show in the 
Table are effective masses obtained from local cosh fits to the correlation function. They are 
labelled by the distance n t at which they are calculated. Although the errors become large 
for n t > 3 (particularly in the confined phase) it is quite clear that there is a discontinuity in 
the value of m s as we cross from the confined to the deconfined phase. (The string operators 
we use have excellent overlaps onto the lightest string states, and we expect ^1] the n t = 2 
effective mass to provide an excellent approximation to the true mass.) Since the transition 
is strongly first order for SU(8), perhaps the real surprise is that this jump (decrease) in the 
spatial string tension is only about 15%. 

From Table IT3*1 we estimate am s in the confined phase to be am s = 1.377(40). Using 
eqn(|41|) with = 2 we then obtain ay/a s (T = T c ) = 0.349(5). We can compare this 
to the 'T = 0' value by interpolating to (3 — 43.965 the values in Table 9 of giving 
a v /cr s (T = 0) = 0.346(2). Thus we see that to a good approximation the spatial string tension 
does not change at all when T increases from to T c - at least for large enough N. 
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We have also calculated am s on the L s = 8 and L s = 6 lattices that were used in Sec- 
tion l6.2l to calculate the T-dependence of the electric gluon mass in the deconfmed phase. The 
extracted loop masses are listed in Table EH In FigJEl we plot the k = 1 spatial loop mass 
against (T/T c ) 2 , since that is the expected behaviour at high T (see eqn (lTD|) ). We plot the 
best linear fits to the data and we clearly observe good evidence for the early onset of the 
quadratic behaviour. We note that the T = intercept corresponds to relative corrections 
that are 0(1 /T 2 ), and which include the string correction. We note some variation of the 
slope with N, but this is close to what we have observed [2] in the variation of T c j^/a with 
TV and so it is presumably an 0(1/N 2 ) correction. 

In FiglHlwe plot the ratios of the /c-loop masses to the fundamental (k = 1) loop masses. 
At higher T, where string corrections become negligible, this just becomes the ratio of the 
correponding spatial fc-string tensions. For comparison we show the predictions of the Casimir 
scaling conjecture ^2] i n eqnfjHJ). (This also arises naturally in a model where the high-T 
plasma contains a gas of adjoint magnetic pseudo-particles It is clear from FigJUJthat 

we have quite good agreement with this conjecture. By contrast, the trigonometric 'MQCD' 
conjecture in eqnfjHJ) is strongly disfavoured. The inclusion of a string correction strengthens 
this conclusion. 



We pointed out in Section 12721 that deconfinement may proceed in several discrete steps, if the 
relevant centre symmetry breaks through a sequence of subgroups. In that case different 
fc-string operators will acquire non-zero expectation values at the corresponding temperatures, 
T c (k). As we remarked in Section E21 one expects, on general grounds, that T c [k) < T c [k = 1). 
Thus the temperature we have identified as T c using the k = 1 Polyakov loop as the order 
parameter, will be the highest of these temperatures. If, then, we calculate the masses of the 
various /c-strings in the confined phase at T = T c (k = 1), any such expectation values will 
become manifest, and we will have direct evidence for the existence of such multiple phase 
transitions. 

We have performed such mass calculations for many of the values of N, lattice sizes and j3 
values discussed in Section |H1 In Table M we summarised some of our 'cleanest' results. These 
calculations were performed without any vacuum subtraction so that a non-zero vacuum 
expectation value would manifest itself in a (nearly) vanishing mass gap. It is clear that none 
of the masses are near zero, demonstrating that in all these case the fc-string operators have 
not developed non-zero expectation values for T < T c {k = 1). This provides our first, quite 
convincing evidence that deconfinement proceeds through a single phase transition. 

Of course the reason the mass calculations in Table 01 are 'clean' is that the spatial volumes 
used are large enough for the confined-deconfined tunnelling to be negligible. One might 
therefore worry that even if there was an earlier phase transition, the corresponding tunnelling 
would be improbable. This is clearly not going to be the case unless the T c (k) are sufficiently 
close, since otherwise the free energy difference between the relevant phases will promote the 
tunnelling probability at T c (k = 1). Nonetheless this uncertainty suggests that we should also 
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look at smaller volumes where there is significant tunnelling. We have indeed done so and 
find no suggestion for non-zero expectation values for any k > 1 string operators. Such a 
study does however suffer from some small ambiguity to do with separating the confined and 
deconfined phases on smaller volumes. 

As an explicit illustration, we consider SU(4). We take the largest lattice, 20 3 5, on which 
we see any tunnelling (within our statistics) and a value of T slightly below T c ; roughly 
T ~ 0.996T C . In FigJTUlwe plot a histogram of the values taken by \lk=i\, the modulus of the 
Polyakov loop averaged over a single lattice field configuration. We see the very large peak 
at low |/| corresponding to the confined phase (I is complex hence the suppression at I — 0) 
and a small peak at larger |/| corresponding to fields in the deconfined phase. If there was a 
separate k — 2 transition then we would expect that some of the fields in the low \lk=i\ peak 
in FigEl would be confined for k = 2 strings and some would be deconfined. To test whether 
this is so we have taken separately all the fields with \lk=i\ < 0.05 and |/jt=i| > 0.05 (i.e. 
mid-way between the confined and deconfined peaks) and for each of these two ensembles we 
plot separately the distribution of \lk=2\ ( m the totally antisymmetric representation). These 
are plotted in Fig^^and we see that there is no evidence at all for a k = 2 deconfined phase 
for which k = 1 strings are still confined. We infer that there is a single deconfining transition 
at any N. 



8 Interface tensions 

We saw in Section 12.31 that at T = T c the interface tension, a c d, determines the rate of 
tunnelling between the confined and deconfined phases through eqn(|12|). In this Section we 
shall use the observed rate of tunnelling to estimate the iV-dependence of a c d, borrowing 
standard methods that have been used in earlier work 2EJ in the case of SU(3). (For a 
different approach, which appeared as this paper was being written, see |2*7j.) We shall start 
by formulating in more detail the expression for tunnelling, with the focus on making explicit 
all the important iV-dependence. 

We use the average Polyakov loop on each lattice field, I, as our order parameter. (In 
practice we will use its modulus, which leads to some minor changes as described in [26J.) Let 
i = c,d label the confined and deconfined phases and let Piil) be the normalised probability 
distribution. If & is the length scale on which the fluctuations of Polyakov loops decorrelate, 
and if V is large enough, then we can expect that making the usual assumption of a single-scale 
Gaussian distribution 



m) = ae d i (42) 

where 

JccJVjf ; coci (43) 

is a harmless approximation for our purposes. The factor of N 2 in eqn ()43|) follows from 
the usual large-iV factorisation arguments for gauge invariant operators, which imply that 
(P) = (Z) 2 (l + 0(l/N 2 )) in both phases. 
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At T = T c we have tunnelling and both Gaussians peaks are present. The probability that a 
value of I mid-way between the peaks comes from one of these Gaussians is oc exp{— b{N 2 V}, 
where the hi are some constants, and this is much less than the tunnelling probability, oc 
exp{— 2a c dA/T}, for large enough V and/or N. So in the latter limit, on an L\L t lattice, the 
observed probability of such a value of / should reflect the tunnelling probability 



p t oc e -2 — t — = e ~ 2a a ° dL « Lt (44) 
in lattice units. 

Let us assume that we have defined T c on a finite volume as the temperature at which the 
probability of the confined phase equals that of the deconfined phase. (The maximum of the 
specific heat or of the loop susceptibility will correspond to this definition at large V.) We 
see from the above that if we consider the overall probability distribution for I then the value 
at the minimum p m i n , between the two Gaussian peaks, is related to the value at one of these 
peaks, p max , by 

Pmax 

The factor of L 2 comes from translations of the two walls (and being careful with the way this 
affects the measure). The fact that there are A" deconfined phases in the Euclidean system 
does not alter this result. Note that we have neglected the contribution from the interiors of 
the domain walls and various factors that have non-zero large- N limits. 

Apart from the obvious caveats, an additional uncertainty in the above analysis is that it 
ignores the contribution of tunnellings between the A" different deconfined phases. If we label 
the deconfined vacuum by k, corresponding to the fc'th element of Z^, then one can see that 
a tunnelling between the k and k' vacua will pass through fields with 

in~|J d |.|cos{7r^^}| . (46) 

For k — k' ~ 0(N/2) the value of |Z| will be partway between \l c \ and \ld\, i.e. precisely where 
we have applied eqn (j43j) . Since we expect (at least at high T where it is calculable [T7j) that 
&kk' ~ 0(N 2 ) when k — k! ~ 0(N/2), such tunnellings are strongly suppressed, and will only 
be important for N, V — > oo if a c d ~ 0(N 2 ) so that the deconfining tunnellings are similarly 
suppressed. In principle one can identify and exclude these deconfined-deconfined tunnellings. 
In practice this requires large enough V for the transitions to be very well separated. This 
represents a potentially significant systematic error. 

As an aside, we remark that the above problem only arises because we have chosen to 
use I as our order parameter. We could use the average plaquette instead since that also 
takes different values in the confined and deconfined phases. However because its value is the 
same in all A^ deconfined vacua, tunnellings between the latter will not populate the range of 
plaquette values that signal tunnelling between confined and deconfined phases. The reason 
we do not analyse the plaquette probability distributions is that the separation between the 
Gaussian peaks is essentially non-existent for our range of V and N. 
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Since our calculations are at T ~ T c it is customary to define the dimensionless surface 
tension 

- W ( 47 ) 

c 

in which case we see from eqn ()45j) that 

^ = ^ + §~A ^ ^ - ^ iV + c} (48) 

Pmax <Z-^ S * 

where c is a constant which may contain large 0(1/N 2 ) corrections. Capillary waves along 
the interface |2Hj produce a constant that can be absorbed into the c in eqn ()48|) . (The next 
correction [2E| is 0(1/L 2 ) and we shall ignore it here.) We now apply the above formalism, 
and calculate a ca - from our calculated Polyakov loop probability distributions. 

In our calculation we choose T c as the temperature at which the peak heights of the confined 
and deconfined Polyakov loop probability distributions are equal. This differs slightly from 
our above discussion but, since both peaks shrink with N and V in the same way, this does 
not affect the argument in a material way. Defining 



J cd 



T 2 n T 

1J t i Fmin ■ Lj 



2L 2 Pmax 4Z/^ 



we see from eqn f!48|) that 



^-§>{c-lnN} + ... 

cr cd . (50) 



In FigH21we plot our calculated values of against (L t /L s ) 2 and show the best fits obtained 
by including the leading large-\^ correction shown in eqn ()50|) . We note that the slope of the 
best fit decreases with increasing N whereas from eqn ()5U|) one would expect the opposite. It 
is clear that higher order corrections in iV or V must be important here. For example, we 
find that the two Gaussian distributions are not well separated for most of our volumes. Thus 
there must be some doubt about the reliability of our laige-V extrapolations. 

Ignoring this uncertainty, we extract the interface tensions from the large-l^ extrapolations 
in Fig^J and list them in Table El Plotting them against N in Fig^] we observe a very 
strong increase with N. The fact that a c d increases by a factor of ~ 30 when we go from 
SU(3) to SU(8) implies, using eqn(|44)l. that we must decrease L s by a factor ~ 5.5 if we want 
to maintain a comparable tunnelling rate when we go from iV = 3 to N = 8. 

To test the functional form of the increase of a c d with N, we show in Fig^Jthe best fits 
corresponding to both linear 

^ = 0.118(3)JV- 0.333(9) ; x 2 /n df = 6.5 (51) 

c 
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and quadratic variations 

= 0.0138(3)iV 2 - 0.104(3) ; X 2 /n df = 2.7. (52) 

c 

While a cursory glance at FiglTHl might have suggested that a linear behaviour looks better, we 
see that this is not supported by the statistical analysis which strongly prefers a c d oc iV 2 . How- 
ever the x 2 is driven by the very accurate smaller- N values, which might be sub-asymptotic. 
So it is clear that more accurate calculations with a better control of the systematic errors are 
needed, particularly for iV > 8, if we are to be able to answer this question. 



9 SU(2) 



The SU(2) gauge theory is linearly confining at low T and deconfined at high T. Therefore 
there will be a phase transition at some T = T c . If the phase transition were first order 
then, in the ^ = 00 thermodynamic limit, dlnZ/dft would be discontinuous at T = T c . This 
is the case, as we have seen, for N > 3. In SU(2), however, it is well known that there 
is no observed discontinuity, and that instead the effective string tension appears to vanish 
smoothly as T — > T~ . The corresponding correlation length diverges as T — > T~ and the 
transition is thus taken to be second order. 

In FigElwe provide an example of how the correlation length varies with /3 and the spatial 
volume on a set of L 3 5 lattices for which V = 00 deconfinement occurs at j3 c = 2.3714(6) (see 
Table 4 of 2\). What is actually shown here is £t = l/am t where m< is the mass of the lightest 
flux loop that winds around the time-torus, and which is related to a e ff(T) by eqn (fT6*j) when 
V — > 00. We see how on a fixed volume near T c the value of £ t is constrained by the spatial 
size, but that it is consistent with diverging at T c as we take V — > 00. This is to be contrasted 
with the values of amt at T = T c for iV > 3 that are listed in Table El (for the same value of 
a in units of T c ). 

The usual way to identify a first or second order transition at (3 = (3 C is to calculate the 
specific heat C(f3) 

P 2 C ^)- V d(3 2 ■ ( 53 ) 

If V is finite but large, then C(j3) should have a peak at some (3 = (3 C (V) v -^? fj c an d this 
peak should diverge as 

C{(3 C {V)) oc V a (54) 

where a = 1 for a first order transition and a < 1 for a second order transition. The behaviour 
for a first order transition follows trivially from the discontinuity in the action, when one notes 
that C is proportional to the integrated connected plaquette correlation function 

C (P) a y\(zZ u pzZ u p'} ~ E m p) 2 oc J2((u p - u)(u P0 -u)). (55) 
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Of course it will only be seen if the Monte Carlo calculation is long enough to sample both 
phases. For a second order transition the divergence of C arises because of the vanishing 
mass, which implies a divergence in the integrated correlator in eqn (j55|) . For the first order 
transition, for N > 3, a growing specific heat peak is precisely what we see. However for the 
SU(2) transition C(/3) not only does not grow with V in the neighbourhood of /3 C , but shows 
no peak at all; despite the fact that, as we see in Fig lT^ there is a rapidly growing correlation 
length for these volumes. This is the implicit reason why earlier studies have not used C(f3) 
to study the SU(2) transition, but have rather used the Polyakov loop susceptibility, which 
shows the expected diverging behaviour. 

We are already familiar with one example of a striking phenomenon in the high T Euclidean 
system which is visible in Polyakov loops but which is not directly reflected in the real thermal 
properties of the gauge theory: the Z N symmetry breaking and vacuum degeneracy discussed 
in Section 12.51 Here we have a diverging correlation length which has only been visible in 
Polyakov loops. Is it possible that this divergence is also not reflected in the thermal properties 
of the gauge theory, so that the phase transition is not second but higher order? This may 
seem very unlikely but provides us with a motivation for understanding better what is going 
on. 

The vanishing of the mass of the timelike flux loop at f3 c implies, from eqn(|16p. the van- 
ishing of the effective string tension, cr e //(T), at T = T c = l/a([3 c )L t . When we calculate the 
overlap of our (unsmeared) Polyakov loop onto this lightest flux loop, at say (3 = 2.3725 ~ (3 C , 
we see that it remains very large, ~ 0.95, on all our volumes. This implies that the effective 
potential between static sources will show a linear behaviour V(r) ~ a e ff(T)r from some 
relatively short distance, and that the potential will flatten from this distance onwards as 
T — > T c . Thus the vanishing of the SU(2) effective string tension at the transition is a robust, 
highly visible feature of the thermal gauge theory. 

Nonetheless this is a property of the gauge theory probed with static fundamental sources, 
and is not related in an obvious way to the vanishing in the pure gauge theory of a screening 
length which couples to the plaquette. Indeed we know that in the confining phase a non- 
contractible winding flux loop must have an exactly zero overlap on to a contractible Wilson 
loop such as the plaquette, and so cannot contribute directly to C((3). But this is not the 
end of the story. Scattering states composed of two such (mutually conjugate) flux loops will 
have energies all the way down to zero at T c and these can have non-zero overlaps onto the 
plaquette. We now estimate this overlap. 

Consider a 20 3 5 lattice at (3 C = 2.3675, which is very close to (3 C . We consider operators 
composed of p = sums of spatial (magnetic) and timelike (electric) plaquettes and also an 
operator that is composed of a product of two p = Polyakov loops. (Since we reserve the 
label t for the short torus, we take correlations in the ^-direction and the momentum p is 
defined in the 3-space orthogonal to the latter.) Looking at the effective mass plot of the 
latter it is clear that at this a this operator has a very good overlap onto the state composed 
of two of the lightest p = timelike flux loops. We will therefore make the approximation 
that this operator is in fact the wavefunctional for this two loop state, \U). We label the 
latter operator by <fiu, and the plaquette sums by u s ,u t respectively. Then the normalised 
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cross-correlations give us the matrix element which determines the desired overlap 



(vac\u a d>u\vac) (vac\u a \ll) 1 0.027(3) : a = s . . 

~ ~ ; (56) 



(vac\u a u a \vac) 2 (vac\(J)ii(f)u\vac) 2 (vac\u a u a \vac) 2 ) 0.042(3) : a = t. 

This is very small, and one naturally wonders if this smallness does not indicate that the 
overlap actually vanishes either as one approaches V = 00 or as one approaches (3 = (3 C . This 
we shall now show is not the case. We have performed calculations on 14 3 5 and 12 3 5 lattices 
which show that the overlap in eqn ()56|) is volume independent, within errors. We have also 
performed calculations at (3 — 2.3600 where T — T c is about 3 times larger than in the above 
example, and again we find that the overlap is unchanged. We thus conclude that at the 
lattice spacing a ~ 1/5T C there is in fact a non-zero overlap between the plaquette and the 
double flux loop whose mass, we recall, vanishes as V — > 00 and T — > T c . That is to say, the 
thermal average of a plaquette-plaquette correlator must, at large enough separations, display 
the contribution of this diverging correlation length, and so the deconfining transition in the 
SU(2) thermal gauge theory is indeed a conventional second order one. 

In the normalised p = plaquette-plaquette correlator, the contribution of the double flux 
loop will be 

C ™< ^ ' 1 " ' _ e -2am t n z + /gjX 



\{vac\ 


u a \ll)\ 2 


{vac\ 


u 




vac) 



C uu (0) 

(ignoring periodicity). We can use this to estimate how large we have to make L s for the 
diverging correlation length to start making a visible contribution. First we note that we can 
replace J2 P u v m e( l n ^^ by the sum over z of the p = sums of plaquettes. And we can make 
a crude estimate for amt(L s ) using FiglTH Using the latter and the overlaps in eqn (J56j) . we 
can easily see, via a very approximate calculation, that for the term in eqn (|57jl to contribute 
about half of the specific heat at /3 = (3 C , we need to have a lattice with L s ~ 2000. It is now 
clear why even on our 'very large' 40 3 5 lattices we saw no sign of a growing specific heat peak! 

The above calculations have been for a specific small range of a close to a = 1/5T C . As 
a — > all physical states decouple from the plaquette, since wavelengths on physical scales 
contribute only ~ cl^Aqcd of its value. So in the continuum one would want to use a suitably 
regularised form of the action when calculating the specific heat. Since we have no reason to 
expect any pathology as a — > we have not attempted a detailed study at smaller a. We have 
however performed some calculations at a = 1/8T C which confirm that the overlap, as defined 
in eqn fl56j) . remains non-zero, and that it is smaller, roughly as one would expect. 

Finally we comment that the above near-decoupling will obviously also occur for N > 2 
where the transition is first order. This means that the specific heat will give a clear, precocious 
signal for the developing two peak structure in the plaquette distribution, unobscured by an 
additional term coming from an increasing correlation length (albeit increasing only to some 
finite value). This is to be contrasted with the Polyakov loop susceptibility. This effect is 
strongest for the weakly first order transition in SU(3) - compare Fig. 5 and Fig. 3 in j2|. 
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10 The bulk phase transition 



It is well known that SU(2) and SU(3) lattice gauge theories with the plaquette action exhibit 
a rapid crossover between the strong and weak coupling regions. This crossover becomes 
more pronounced for SU(4) [2HJ and, as we have found in the present work, transforms into a 
weakly first order transition for SU(5) and becomes robustly first order for SU(iV > 6). The 
strong coupling bulk phase is believed to contain vortices and monopoles which disorder 
Wilson loops down to the ultraviolet length scale, so that the string tension in lattice units 
is a 2 a ~ 0(1). On the weak coupling side, the physics in the ultraviolet is determined by 
asymptotic freedom and so a 2 a ~ 0(a 2 ) <C 1. Only calculations on the weak coupling side 
are useful for the continuum extrapolation. 

The N — > oo limit has been studied within the twisted Eguchi-Kawai formalism, and the 
critical inverse bare 't Hooft coupling is [SD] 

X^(N = oo) = Jim = 0.3596(2). (58) 

In Table we list the values of (3b thai we have obtained for 5 < iV < 12 during the course 
of our calculations. For iV = 5 we can estimate (3b directly but for N > Q the interface tension 
becomes large enough that all we see is a hysteresis effect. Once the interface tension is large 
the size of the hysteresis varies weakly with both the volume size, as we see in Table and 
also with the length of Monte Carlo run at each value of (3 (within practical limits). This 
leads to a large uncertainty in (3 b , which can be eliminated by adding a suitable function of 
the order parameter to the action, as in jHOj- However our purpose here is not to study the 
bulk transition per se but rather to determine its influence, if any, on our calculations of T c . 
Nonetheless we note that our results for N > 6 are consistent with eqn (|58j) together with 
a modest 0(1/N 2 ) correction. (Note that this will not be the same 0(1 /N 2 ) correction as 
determined in [3U] since the Eguchi-Kawai equivalence only exists at the leading planar level.) 

We compare the values of (3b in Table HH1 to the values of (3 C in Tables HE We note 
that for SU(6) the deconfining transitions for L t > 5 are safely beyond the bulk transition. 
Interpolating to SU(5) we would expect the L t = 5 deconfining transition to be at (3 c (L t = 
5) ~ 16.87, which is again well beyond (3b- For SU(8) on the other hand, the hysteresis has 
widened and now encompasses both (3 c (L t = 5) and (3 c (L t = 6). These calculations of (3 C have 
been performed by working on the weak coupling branch of the hysteresis curve and clearly 
we need to know if this introduces a bias. 

A useful order parameter for the bulk transition is the average plaquette u p . For larger 
N its discontinuity is Au p ~ 0.13. Thus, in contrast to the deconfining transition, where the 
discontinuity is only 0(a 4 ), here it is 0(1) and so even partial tunnellings are easily visible 
on single lattice fields. We have examined the plaquette distributions at (3 ~ (3 c (L t = 5) and 
we find no trace of even a partial attempt at tunnelling. It is clear that the barrier between 
the bulk and weak coupling phases is so large that the deconfining phase transition in the 
weak coupling phase at a ~ 1/5T C is largely unaffected by the fact that it occurs in only a 
local rather than in the absolute minimum of the effective potential - except perhaps for some 
minor distortion of the potential that we have earlier hypothesised might be the reason for 
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the unexpectedly large lattice spacing corrections to the latent heat at this a. 

11 Conclusions 

In this paper we studied various properties of SU(iV) gauge theories at and around the de- 
confining phase transition. We checked that there is indeed only one transition, i.e. that the 
centre symmetry is not broken in steps, with each step corresponding to the condensation of 
an appropriate /c-string. We showed that the transition is robustly first order at large N with 
a latent heat Lh oc iV 2 as expected, and an interface tension a c d that grows with N, although 
our calculations were unable to determine whether this growth is oc iV or oc iV 2 . In the former 
case the transition would become a real phase transition even on a small volume as iV — > oo 
(which is possible because the number of degrees of freedom per unit volume diverges in that 
limit). In the latter case one expects a growing hysteresis, even on small volumes. If this 
hysteresis is large enough one might encounter the 'Hagedorn' string condensation transition 
along the confining branch, above T c . 

We calculated the spatial string tension and showed that at larger iV its value in the 
confining phase appears to be entirely independent of T. However, when we compare the 
spatial string tensions in the two phases at T = T c it is clear that there is a significant 
discontinuity. This is easy to see for larger N because the transition is robustly first order 
and the growing interface tension ensures that even on moderately large lattices there is no 
tunnelling, or even attempted tunnelling. So one can study the properties of the two phases 
separately at T = T c . For example, on a 12 3 5 lattice for SU(8). By contrast, in the well-studied 
case of SU(3), even on a 64 3 5 lattices we still observe tunnelling. Thus such a discontinuity 
has not been observed in SU(3) although it undoubtedly exists there as well. 

At higher T we have seen that the spatial string tension grows oc T 2 and that it does 
so essentially from T = T c for larger N. The spatial fc-string tensions are consistent with 
Casimir Scaling. For high T the Debye mass grows oc T as expected, although there are large 
deviations near T c ,. This is not surprising because the leading-order theoretical expression, 
m 2 D oc g 2 (T)T 2 contains an additional 'logarithmic' dependence on T in the coupling. 

We also showed why the SU(2) transition has no peak in the specific heat, let alone one 
that grows with V, although there is a diverging correlation length. We explicitly showed that 
this correlation length will eventually appear as a visible screening length in thermal averages, 
so that the transition is a bona fide second order one. For larger N we discussed the bulk 
transition, in the shadow of whose hysteresis we have performed some of our deconfinement 
studies. 

The main motivation for studying SU(iV) gauge theories at large iV is the hope that 
SU(3) might be 'close' to this limit and the gauge theory is expected to be much simpler at 
N = oo than at iV = 3. The former has been largely confirmed by recent lattice calculations 
[H3 El 12] ■ Signs of the latter are found throughout our calculations in this paper; for example 
the reduction in finite volume effects at the deconfinement transition, which naturally points [2 
to a form of Eguchi-Kawai space-time reduction |3T]. On the other hand, while the suppression 
of 'small' fluctuations as N — > oo does indeed lead to the simple idea of a Master Field |3^j . 
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what we have also seen is that there is actually a plethora of Master Fields - confining, 
deconfining, bulk, one for each of the N degenerate phases obtained if a spatial torus is 
reduced below a critical size, not to mention the N 6*-vacua branches jHSl and conjectured 
non-analyticities in scale size [B3]. All this [23J serves to emphasise how much there is still to 
be understood about the physics of this theory. 
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SU(4) 



L t 


L s 


h 


Pc(V = oo) 


T c /y/a 


5 


12-20 


0.090(17) 


10.63727(53) 


0.6148(13) 


6 


16 


0.100(22) 


10.7898(16) 


0.6166(24) 


8 


24 


0.123(27) 


11.0880(22) 


0.6310(30) 



Table 1: Critical values of (3 calculated on L\L t lattices and extrapolated to V = oo using 
eqn(|26|) and the values of h shown. The corresponding values of T c j ' ■sja are shown. 



SU(6) 



L t 


L s 


h 


Pc(V = oo) 


Tjy/a 


5 


8-14 


0.112(19) 


24.5139(24) 


0.5894(36) 


6 


16 


0.128(29) 


24.8467(30) 


0.5956(29) 


8 


16 


0.168(39) 


25.4782(64) 


0.6024(29) 



Table 2: Critical values of (3 calculated on L\L t lattices and extrapolated to V = oo using 
eqn()26|) and the values of h shown. The corresponding values of T c /^/a are shown. 



SU(8) 


L t 


L s 


h 


/3 C (V = oo) 




5 


8,10 


0.134(73) 


43.982(14) 


0.5819(41) 


6 


8 


0.155(84) 


44.535(37) 


0.5850(70) 


8 


10,12 


0.218(123) 


45.654(32) 


0.5888(69) 



Table 3: Critical values of (3 calculated on L\L t lattices and extrapolated to V = oo using 
eqn(J2H} and the values of h shown. The corresponding values of T c /^/a are shown. 



SU(iV) : continuum 


N 




X 2 /n d f 


2 


0.7091(36) 


0.28 


3 


0.6462(30) 


0.05 


4 


0.6344(42) [(81)] 


3.7[1.0] 


6 


0.6101(51) 


0.02 


8 


0.5928(107) 


0.003 



Table 4: T c /y/a extrapolated to the continuum limit for various SU(iV) gauge theories, with 
the x 2 P er degree of freedom of the best fits. The SU(2) and SU(3) values are taken from 
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A T 

JM 


P > 


P < 


Po 


c 


Cl 


c 2 


c 3 


2 / 


2 


2.1768 


2.5115 


2.3726 


0.2877 


0.9195 


-3.2961 




0.55 


3 


5.6925 


6.3380 


5.8945 


0.2607 


2.6312 


12.2494 


24.0 


0.54 


3 


5.6925 


6.3380 


5.8945 


0.2610 


2.7271 


10.6846 




1.1 


4 


10.637 


11.400 


10.789 


0.2702 


3.5393 


18.7900 




0.9 


6 


24.500 


25.452 


24.670 


0.3082 


4.9649 


34.5288 




0.5 


8 


43.85 


45.700 


44.350 


0.3022 


5.0812 


42.0732 




1.2 



Table 5: Fit parameters to the string tension as a function of (3 using the interpolation 
function in eqn(|33j). 



N 


u 


L s G 


A(u p ) 


3 


5 


20-40 


0.00084(2) 


4 


5 


14-20 


0.00198(5) 


6 


5 


8-14 


0.00238(6) 


8 


5 


8-10 


0.00251(34) 



Table 6: The jump in the average plaquette, A{u p ), between the confined and deconfined 
phases at (3 c {a = 1/5T C ), as obtained from a V — > oo extrapolation of the specific heat peak 
from the range of lattice sizes shown. 



L h at a = l/L t T c 



N 


u 


L s 


P 


A{u p ) 


d In a /dp 


L h /T* 


3 


5 


64 


5.799 


0.00078(3) 


-2.070(9) 


1.413(55) 


4 


5 


32 


10.635 


0.00187(4) 


-1.2878(90) 


5.39(13) 




6 


16 


10.780 


0.00064(5) 


-1.1666(80) 


4.06(34) 




8 


24 


11.085 


0.000165(20) 


-0.9400(90) 


4.30(53) 


6 


5 


16 


24.515 


0.00248(3) 


-0.6347(120) 


14.69(35) 




6 


16 


24.845 


0.00084(2) 


-0.5472(80) 


11.90(35) 




8 


16 


25.46 


0.000202(7) 


-0.4012(76) 


12.00(49) 


8 


5 


12 


43.965 


0.00261(4) 


-0.3717(54) 


25.67(65) 




6 


8 


44.45 


0.00075(5) 


-0.3240(26) 


16.0(1.4) 




8 


12 


45.50 


0.000195(10) 


-0.2304(48) 


17.8(1.1) 



Table 7: The jump in the average plaquette, A(u p ), and the latent heat, L h , at T = T c , for 
the gauge groups and parameters shown. 
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L 


h : continuum 


N 


Lte 


Ll/Tc 


4 


6,8 


1.47(10) 


6 


6,8 


1.87(5) 


8 


6,8 


2.12(9) 



Table 8: The latent heat in physical units after a continuum extrapolation from the L t = 6, 8 
values. 







arrik 


at a ~ 1 /5T C 


: confined 




N 


L s 


P 


k = 1 


k = 2 


k = 3 


3 


64 


5.799 


0.105(8) 






4 


32 


10.635 


0.1967(89) 


0.335(13) 




6 


16 


24.515 


0.2343(87) 


0.474(25) 


0.657(36) 


8 


12 


43.965 


0.2835(103) 


0.485(43) 


0.69(10) 



Table 9: The masses of the lightest time-like /c-loops in the confined phase at T = T c . 







arrik at a ~ 1/5T C 


: deconfined 


N 


L s 




fc = 1 


k = 2 


fc = 3 


3 


64 


5.799 


0.080(6) 






4 


32 


10.635 


0.191(34) 


0.185(32) 




6 


16 


24.515 


0.359(14) 


0.346(13) 


0.346(12) 


8 


12 


43.965 


0.418(14) 


0.405(13) 


0.402(12) 



Table 10: The lightest non-zero masses coupling to the time-like fc-loop in the deconfined 
phase at T = T c . 



am t at a ~ 1/5T C : SU(4) 


L s 


P 


C:k = l 


C: k = 2 


D 


32 


10.635 


0.1967(89) 


0.335(13) 


0.185(32) 


20 


10.635 


0.203(10) 


0.312(20) 


0.215(20) 


16 


10.635 


0.195(20) 


0.24(4) 


0.25(3) 


20 


10.642 






0.240(20) 


20 


10.633 


0.185(10) 


0.255(20) 




16 


10.630 


0.204(15) 


0.27(4) 





Table 11: The lightest non-zero masses coupling to the time-like fc-loop in the confined (C) 
and deconfined (D) phases at T ~ T c on various volumes. 
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2am D for T >T C 


N 




T = T C 


T = 5T c /4 


T = 5T c /3 


T = 5T c /2 


T = 5T c /2; L s = 6 


3 


5.80 




0.615(21) 


1.040(15) 


1.596(40) 




4 


10.635 




0.660(8) 


1.053(19) 


1.643(40) 


1.578(30) 


8 


44.0 


0.527(11) 


0.670(12) 


1.054(20) 


1.599(46) 


1.62(5) 



Table 12: The lightest non-zero mass that couples to the time-like flux loop in the deconfmed 
phase for a ~ 1/5T C on L s = 8 lattices (unless otherwise indicated). 



SU(8) : am s (t = an t ) 


n t 


D 


C 


1 


1.326(7) 


1.527(11) 


2 


1.187(24) 


1.377(40) 


3 


1.173(81) 


1.16(15) 


4 


1.44(39) 


1.20(42) 



Table 13: Effective masses extracted at t — an t from correlators of spacelike Polyakov loops 
on 12 3 5 lattices at T = T c in the confined, C, and deconfmed, D, phases. 



am s (k) for T > T c 


N 


k 


T = T C 


T = 5T c /4 


T = 5T c /3 


T = 5T c /2 


T = 5T c /2; L s = 6 


3 


1 




0.753(7) 


1.119(11) 


2.027(42) 




4 


1 




0.835(7) 


1.192(12) 


2.211(66) 


1.607(23) 




2 




1.126(17) 


1.59(4) 


2.52(32) 


2.09(8) 


8 


1 


0.695(9) 


0.910(11) 


1.311(19) 


2.39(13) 


1.805(32) 




2 


1.171(26) 


1.610(34) 


2.12(12) 




3.1(7) 




3 


1.507(41) 


1.933(81) 


2.29(32) 








4 


1.603(72) 


2.41(23) 


2.75(99) 







Table 14: The spacelike A;-loop masses on 8 3 5 lattices (unless otherwise indicated) in the 
deconfmed phase. 
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interface tension 



N 


<WT C 3 


fit range 


X 2 /n df 


n df 


3 


0.0200(6) 


16-40 


1.51 


3 


4 


0.1208(56) 


14-20 


0.37 


2 


6 


0.394(11) 


8-14 


6.02 


2 


8 


0.56(10) 


8-10 








Table 15: The tension of the confining-deconfining interface, in units of T c , for various SU(iV) 
gauge theories, extracted from the V — > oo extrapolations in Fig^J All for fixed a: aT c = 1/5. 
The best fit \ 2 an d the number of degrees of freedom, njf, are shown. 



bulk transition 


N 


lattice 






5 


8 4 


16.6565(5) 


16.6565(5) 


5 


10 4 


16.6550(5) 


16.6550(5) 


5 


12 4 


16.6550(10) 


16.6550(10) 


6 


4 4 


24.100(25) 


24.475(25) 


6 


6 4 


24.290(5) 


24.480(5) 


8 


4 4 


43.25(5) 


44.85(5) 


8 


6 4 


43.60(5) 


44.90(5) 


12 


4 4 


98.5(1) 


103.9(1) 


12 


8 3 5 


98.7(1) 





Table 16: The location of the strong- weak coupling bulk transition when decreasing or in- 
creasing f3 from large or small values respectively. For various SU(AQ groups on the lattices 
shown. 
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Figure 1: The value of the 't Hooft coupling on the scale a, as obtained from /3 in eqn ()25|) . for 
N = 2(o), 3(D), 4(*), 6(+), 8(«), plotted against the values of a expressed in physical units. 
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Figure 2: The SU(iV) continuum deconfining temperature in units of the string tension, with 
an extrapolation to N = oo using a leading 0(l/N 2 ) correction. 
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Figure 3: The latent heat calculated at a = 1/5T C for various SU(iV) groups. 
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Figure 4: The deconfinirig latent heat in SU(4), o, SU(6), □ and SU(8), o gauge theories at 
various a, with extrapolations to the continuum limit using a 0(a 2 ) correction. 
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Figure 5: The masses of the lightest k — 1 (•) and k — 2 (o) flux loops that wind around the 
time-torus, in the confining phase at T = T c , a = 1/5T C for various SU(iV) gauge theories. 
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Figure 6: The lightest mass that couples to the vacuum subtracted Wilson line winding around 
the time-torus, in the deconfined phase at T = T c and a = 1/5T C for various SU(iV) gauge 
theories. 
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Figure 7: The electric screening mass plotted against T/T c for SU(3) (o) and SU(8) (•). 
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Figure 8: The mass of the lightest k — 1 flux loop that winds around a spatial torus. Plotted 
against T 2 for SU(3) (*), SU(4) (o) and SU(8) (•), all at a ~ 1/5T C . 
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Figure 9: The mass ratio of the k = 2 to k = 1 spatial loops in SU(4) (•) and in SU(8) (o), 
and of the k — 3 to k — 1 loop in SU(8) (•*■). All in the deconfined phase and for a ~ 1/5T C . 
High T Casimir scaling predictions are shown for comparison. 
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Figure 10: Histogram of timelike k — 1 Polyakov loops averaged over single 20 3 5 lattice fields 
in SU(4) at T ~ T c . Separate confined and deconfined peaks are visible. 
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Figure 11: As in FiglTUlbut for (antisymmetric) k = 2. From lattice fields with |Z&=i| < 0.05, 
• , and with \h=i\ > 0.05, o. 
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Figure 12: The quantity S C( ^ defined in eqn fl50|) whose L s = oo intercept is the interface 
tension, a cd /T c 3 . Calculations for SU(3), □, SU(4), O, SU(6), A and SU(8), o. All are for 
a ~ 1/5T C . Linear extrapolations to infinite V are shown. 
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Figure 14: The correlation length near the SU(2) deconfinement transition, as a function of 
the spatial lattice size. For (3 = 2.360 (•),/? = 2.3675 (o) , and (3 = 2.3725 (*) on L^5 lattices 
where p c (V = oo) = 2.3714(6) 0. 
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